Electromagnetic response of broken-symmetry nano-scale clusters 
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A microscopic, non-local response theory is developed to model the interaction of electromagnetic 
, radiation with inhomogeneous nano-scale clusters. The breakdown of classical continuum- field Mie 

theory is demonstrated at a critical coarse- graining threshold, below which macroscopic plasmon 
resonances are replaced by molecular excitations with suppressed spectral intensity. 

The prevalent classical model describing the interaction of visible and infrared electromagnetic radiation with 
nano-scale metallic clusters is based on Mie theory 1]. This local continuum field model which employs empirical 
values of a bulk material's linear optical response has been used to describe plasmon resonances in nanoparticles 
USB However, such a semi-empirical continuum description necessarily breaks down beyond a certain level of 
coarseness introduced by atomic length scales. Thus, it cannot be used to describe the interface between quantum and 
macroscopic regimes. Moreover, extensions of Mie theory to inhomogeneous cluster shapes are commonly restricted to 
low order harmonic expansions (e.g. elliptical distortions) and so do not exhaust the full realm of possible structures. 
In addition, near-field applications, such as surface enhanced Raman scattering 5], are most naturally described 
using a real-space theory that includes the non-local electronic response of inhomogeneous structures, again beyond 
| the scope of Mie theory. 

Ignoring the limitations of Mie theory can lead to conclusions of limited validity and usefulness. For example, it 
was recently argued 4] that resonances in the long wavelength limit are exclusively controlled by the object's shape, 
but otherwise scale-invariant, so that scaling of nanoparticles would not alter their characteristic plasmon frequencies. 
In particular, it was stated that an infinitely long elliptic rod, whose surface in Cartesian coordinates is defined by 
(x/a) 2 + (y/b) 2 = R 2 , has two geometric scalable plasmon resonances that depend only on the ratio between semi-axis 
lengths a and b. In this work we develop a microscopic model that demonstrates the breakdown of this concept at 
atomic scales, whereas for large cluster sizes the classical predictions for the plasmon resonances are reproduced. We 
then show how this approach is useful in describing the coexistence of excitations of quantum and classical character 
f^. . in inhomogeneous structures on the nanoscale. 

To capture the main single-particle and collective aspects of light-matter interaction in an inhomogeneous nanoscale 
system we adopt the well established linear response approximation 6]. Starting from the Schrodinger equation for 
noninteracting electrons with mass m e and charge e moving in potential V(r) we have 



H^i(r) = { -^—V 2 + V(r) ) *<(r) = M(r). (1) 
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Eq. [T] is solved simultaneously with the Poisson equation that determines the local potential due to the spatial 
distribution of the positive background charges. Using the jellium approximation, the resulting potential is implicitly 
given by \7 2 V(r) = 4e7rp(r), where the density of the positive background charge p(r) satisfies the condition of 
neutrality inside the nanostructure so that J p(r)dr = iV e i, where N e \ is the number of electrons. 

Within linear response theory 0, El ? the induced electron charge density due to an external field of frequency uo is 
given by p[ n d( T ^) = J x( r ? Y ' •> ^O^totC 1 "'? oS)dr' ^ where the total self-consistent potential is determined by 0t o t( r> — 
0ext( r/ >^) + ( / ) ind( r/ 7 ^) • Because the significance of non-local effects for inhomogeneous systems, including spatial 
dispersion of the dielectric function, is widely recognized 0,0? f° r an inhomogeneous cluster of arbitrary shape we 
use the non-local density-density response function x(r, r', uo) that may be determined from a real-space representation 
0. Within the random phase approximation 8] one finds x(r,r',u;) = E E,-E~-hu-\ ^*( r )^( r 0^K r ')^j( r )' 

where f(Ei) is the Fermi filling factor and the small constant 7 describes level broadening. Electron eigenenergies 
Ei and eigenfunctions are obtained numerically using a 6th order real-space discretization of the Schrodinger 
equation in three dimensions with a star-like stencil on a simple cubic lattice. Unless explicitly stated otherwise, zero 
Derichlet boundary conditions are used to calculate the electron wave functions in this paper. The large scale sparse 
eigenproblem is solved using ARPACK . The induced potential 0i n d is determined from the self-consistent integral 
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equation 

s- d(r a,) - v m)-m) p?w 

c^r, lo) Ei-Ej -tiw-ij J |r - r'| 

x | **(r'> tot (r>)^(r")dr", (2) 

and the induced field is found via E in d(r, uj) = — V(/>i n d( r , u). For sufficiently small structures, when the typical 
wavelength of the external radiation is larger than the structure, retardation effects can be ignored. Note, that the 
induced field is not necessarily collinear with the applied external field E ext (r). 

The integral equation Eq.© is discretized on a real-space mesh, using a simple cubic lattice with lattice constant 
L. The natural energy scale Eq is defined by Eo = h 2 /(2m e L 2 ). For example, in the case of GaAs ID nanorods [l2| 
with conduction band carrier density 10 18 cm -3 the discretization lattice length scale is L = 1.28 nm, the density 
p = 1/480 in units of L~ 3 and the energy scale Eq = 330 meV. 
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FIG. 1: Logarithm of the energy of the induced electric field (log 10 (Wi n d)) ( see text) in an elliptic metallic rod due to applied external 
field Eext along the (1,1,0) direction as a function of face surface size R and photon energy Tvjj of the external field in units of Eq. Aspect 
ratio (semi-minor to semi-major axis) a : b = 1 : 1.3. Cross section of the elliptical rod and the vector of the applied external field E ex t are 
depicted in the lower right corner. Temperature T = K, carrier density p = (1/480) L~ 3 , and 7 = 10 -3 E . All real space calculations 
are performed on a mesh of simple cubic geometry with lattice constant L. 

We illustrate the application of this approach to an infinitely long elliptic rod illuminated by an external field. The 
rod is aligned in the z direction and the electric field polarization is along the (1,1,0) direction. For the wave functions 
we assume periodic boundary conditions along the z direction. To quantify the response of arbitrary geometries to 
the applied field, we calculate the energy of the induced field, defined by W ind (u) = (1/2) / |E ind (r,^)| 2 dr. In Fig. □ 
logarithm of the energy of the induced field (log 10 (Wmd)) is displayed as a function of the applied external field photon 
energy hw and the characteristic system size R for the aspect ratio a : b = 1 : 1.3. For sufficiently large rod sizes, one 
clearly observes two plasmon resonances, cj + and cj_, consistent with earlier predictions based on Mie theory [j]. Our 
method confirms that these resonances occur at u+ = co p (b/(a + b)) 1 / 2 and uj- = co p (a/(a + 6)) 1 / 2 , where co p is the 
bulk plasmon frequency uj p — y^^ire 2 p/m e . The grey scale in Fig. is on a logarithmic scale, demonstrating that the 
spectral intensity for large rod sizes is orders of magnitude greater compared to smaller rod sizes. Most importantly, it 
is evident from this figure that the classical picture of two well-defined resonances breaks down below a characteristic 
system size. For sufficiently small rod sizes, the two macroscopic resonances split into multi- level molecular excitations, 
with the overall spectral weight shifting towards lower energies. For the chosen parameters this transition occurs at 
R c ~ 6.5 L. Below R C1 the spectral intensity of the energy levels is reduced because fewer electrons participate in 
the individual resonances. The physics determining the value of R c may be illustrated considering an infinitely long 
cylindrical rod of radius R and electron density p. In the classical regime, the observed collective oscillations are only 
weakly damped, indicating that they are well separated from the quasi-continuum of single-particle excitations. This 
leads to the condition uo^R/vy ^> 1, stating that the plasmon phase velocity is greater than the Fermi velocity vy of 
the electrons For r < R the electrons are trapped in a harmonic potential due to the uniform positive background 
|l4j and the characteristic collective frequency is uo p = ^/A-Ke 2 p/ (2m e ). Estimating the Fermi velocity using a bulk 
value v F = (37r 2 p) 1 / 3 fi/m e one obtains R^> R c = tt 17 ^ 173 — _^ _ ~ 3 4 £ for an electron density p = (1/480) L -3 , 
which is in reasonable, but only approximate, agreement with the calculated threshold R c w 6.5 L shown in Fig. 1. 
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FIG. 2: (a) Logarithm of the energy of the induced electric field (log 10 (Wi n d)) m an elliptic rod of size R = 11 L with aspect ratio 
a : b = 1 : 1.3 as a function of photon energy fru; of the external field in units of Eq. Carrier density p = (1/480) L~ 3 and 7 = 10 -3 
Eq. The direction of the external field is indicated by the arrow. Inset: induced charge density at the resonant frequencies o;+ and uj—. 
The boundary of the rod is shown using the solid line, and the dotted line shows the set of classical turning points, corresponding to the 
positive background potential, (b) Same as (a), but for size R = 3 L. 



For system sizes below R c the dominant excitations are observed to shift towards lower energies as the positive 
background potential becomes increasingly anharmonic. The harmonicity criterion for excitations can be expressed 
as (Jq/R < 1, where <Jo = ) 1//2 is the classical turning point for an electron in the ground state of the harmonic 

potential. This is equivalent to the condition R> R c = ( 2 m pe 2 yielding R > 3 L. 

The nature of the excitations changes as the characteristic size R crosses the quantum threshold R c . In Fig. we 
show logarithm of the calculated energy of the induced field (log 10 (Wmd)) as a function of external field frequency. As 
shown in Fig. 2(a), for relatively large rod sizes, the two distinct plasmon resonances labelled ou- and uj+ correspond 
to two orthogonal bipolar charge distributions, as indicated in the inset. The spatial orientations of these induced 
resonances are aligned with the semi-major and semi-minor axes, and do not depend on the direction of the incident 
field. In contrast, for system sizes less than R c (Fig. GJb)), the excitation spectrum consists of several lower-intensity 
modes, dominated by a low frequency resonance. However, the spatial characteristics of these modes are different 
from the classical limit, i.e. they are aligned with the incident field, indicating that Mie theory breaks down in the 
quantum regime. 

An advantage of the self-consistent non-local response theory described in this work is that it extends naturally to 
inhomogeneous structures with nano-scale and atomically sharp edges and corners. This is of great practical interest 
to nano-photonic applications since such features are commonly associated with substantial local field enhancements 
PJ. However, the intuition driving such expectations is usually derived from classical continuum field theory which 
may not be applicable in this regime. When investigating structures at the microscopic level one needs to account 
quantitatively for changes in the response spectrum. One would also like to describe the breakdown of the continuum 
picture and quantum-mechanical discretization effects which may ultimately lead to new and interesting functionalities 
that are accessible by nano-scale engineering. 

To explore this we consider a representative asymmetric nanostructured system of material volume V ma t = 5516 
L 3 consisting of a conjoined sphere and rectangular bar (Fig. EH a)) whose response strongly depends on the carrier 
concentration (Fig. E^b)). Plots of induced charge distributions in response to an external electric field are shown for 
different carrier concentrations and frequencies of the incident radiation. 

For N e \ = 100 in the low- frequency limit (Fig. Efc)), this leads to a dipole-like anisotropic response for which the 
induced field is not collinear to the external field E ext (indicated by the arrow). Note, that in this case the induced 
charge is localized within the boundaries of the nanostructure, and follows the nanostructure's physical bounds given 
by the positive charge distribution, in agreement with expectations from classical field theory. At this relatively high 
carrier concentration and at a high external field frequency of uj = 0.415 Eq, one observes a complex local enhancement 
of the induced charge density (Fig. Efd)). For the low carrier concentration (N e \ = 5), shown in Fig. OJe) and (f) the 
induced charge density arises from excitation of low energy eigenstates and the overall response is weaker. Note, the 
induced charge is not localized within the boundaries of the nanostructure. Varying the carrier concentration hence 
acts as a "switch" that can activate different resonant regions in a broken-symmetry atomic-scale structure. 

At higher frequencies, the system response can be even more complex. For example, in the case of low carrier 
concentrations and resonance frequency ou = 0.085 Eq, the induced dipole field is rotated with respect to the static 
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FIG. 3: (a) Illustration of three-dimensional conjoined sphere and rectangular bar. (b) Logarithm of the induced field energy in the 
structure for iV e i = 100 (solid line) and iV e i = 5 (dotted line) as a function of the external field frequency. The material volume is 
V^nat = 5516 L 3 . The arrows in the figure show the frequencies chosen for results of calculations shown in (d) and (f). Induced charge 
density for a nanoscale asymmetric structure for the indicated charge densities and frequencies of the external field. Temperature T = 
K and 7 = 10 -3 Eq. (c) N e \ = 100, huj = 0, the direction of the external field is indicated by the arrow; (d) N e \ = 100, faj = 0.415 Eq. 
(e) N e i = 5, fa = 0. (f) N eY = 5, TvjJ = 0.085 E . 




FIG. 4: (a) Three-dimensional image of a structure consisting of a sharp tip and two parallel flat plates of material volume Vmat = 33404 
L 3 containing 7V e i = 207 electrons. The incident electric field is along the (0, 0, 1) direction, temperature T = K and 7 = 10 -3 Eq. (b) 
Logarithm of the induced electric field energy in the structure as a function of the field frequency, the arrow shows the frequency chosen for 
results of calculations shown in (d). (c) Logarithm of the intensity of the induced electric field (log 10 (|E i n d( r )| 2 )) m static limit, fru; = 0. 
(d) At resonance, ftu = 0.135 i^o- 



limit (Fig. EJf))- Hence, by tuning the frequency of the external field, discrete quantum states within anisotropic 
nanostructures can either be accessed or avoided. This control exposes various quantum functionalities that are 
beyond conventional Mie theory. 

As an example, let us consider the induced electric field at the interface between a sharp tip and a flat surface 
illuminated by plane- wave radiation incident in the z direction (Fig. Ufa)). The total material volume of the system 
is V mat = 33404 L 3 . 
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In the static limit, a moderate enhancement of field intensity is found in close proximity to the tip. Results of 
calculating logarithm of the induced field intensity (log 10 (|E ind( r )| 2 )) are shown in Fig. Etc). While most of the 
induced charge density accumulates at the surface, there is considerable penetration into the bulk. In contrast, when 
the same structure is illuminated at resonance, the induced plasmonic response is much more intense than the low 
frequency limit. Here, the induced field intensity increases by 2 orders of magnitude, while the spatial dispersion 
of the "hot spot", as measured by full width half maximum, remains approximately the same (see Fig. Efd)). 
This demonstrates that custom-designed atomic-scale tip structures may be used to control the near field and fully 
quantum-mechanical modeling can quantitatively account for local enhancements, as well as screening and focusing 
effects. 

In summary, a non-local response theory describing the interaction of electromagnetic radiation with inhomogeneous 
nano-scale structures has been developed and implemented. This approach may be used to describe metals and 
semiconductors in the metallic regime. For simple geometries semi-empirical continuum-field Mie theory is found 
to break down beyond a critical coarseness. At this level of coarse graining, the response is no longer exclusively 
determined by particle shape. In more complex inhomogeneous geometries excitations of quantum and classical 
character can coexist. For objects with nano-scale sharp features, our real-space response theory uncovers new 
functionalities, such as local resonances that are activated by tuning the carrier concentration or the frequency of the 
incident field. 

We acknowledge discussions with P. B. Littlewood and V. Kresin. This work is supported by DARPA. Computa- 
tional facilities have been generously provided by the HPC Center at USC. 
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